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We address a problem of generating a robust entangling gate between electronic 
and nuclear spins in the system of a single nitrogen-vacany centre coupled to a near- 
est ^'^ Carbon atom in diamond against certain types of systematic errors such as 
pulse-length and off-resonance errors. We analyse the robustness of various control 



^ ' schemes: sequential pulses, composite pulses and numerically-optimised pulses. We 

(N : 

find that numerically-optimised pulses, produced by the gradient ascent pulse en- 
1^ ' gineering algorithm (GRAPE), are more robust than the composite pulses and the 

sequential pulses. The optimised pulses can also be implemented in a faster time 
^ ■ than the composite pulses. 
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devices. Hence, theoretical and experimental investigations into the control of quantum in- 
formation in such systems have progressed quite rapidly and candidate technologies include 
r'hosphorus donor electrons in Silicon nanostructures 1, ^, Gallium Arsenide quantum dots 
3|, superconducting single-Cooper-pair boxes [J], circuit quantum electrodynamics [5], and 
single nitrogen- vacancy (NV) centres in diamond jsl. A system based on NV centres in 
diamond is attractive because at room temperature, the NV centre displays strong spin po- 
larisation under optical pumping, exhibits remarkable photostability and shows high fluores- 
cence quantum yield Q]. The electronic spins of the NV centre are initialised and measured 
by optical means and can be manipulated by a microwave radiation. Through a hyperfine 
coupling, an interaction between the electronic spins of the NV centre and a nuclear spin 
of a nearby ^'^C atom can be exploited to encode two-qubits and quantum logic can be exe- 



cuted via the application on microwave and radio-frequency radiation. Oservations of Rabi 
oscillations of the single electronic spin j^, and the single nuclear spin 9] in the NV centre, 
performed by optically detected magnetic resonance (ODMR) techniques, have paved the 
way for a realisation of the NV centre based quantum computer A more advanced 

quantum information processing (QIP) task, that is a demonstration of three-qubit entan- 
glement has been recently reported in the literature llj. This latest development points 
to possibilities of multi-qubit QIP implementations. However, the use of accidental nearby 
nuclear spins poses some difficulties towards scaling up this technology. 

Because the qubit system interacts with an environment, QIP is subject to unavoidable 
errors. Hence, error control and avoidance schemes are necessary to achieve reliable quantum 
computation. The errors can occur either in a random or systematic fashion. The random 
errors are due to decoherence processes while systematic errors occur when the physical 
apparatuses controlling the dynamics of the system operate in an imprecise but reproducible 
manner. In the spin system driven by microwave or radio-frequency radiation, non ideal 
values of the radiation properties (i.e. amplitude, phase, duration and frequency), cause such 
systematic errors. A common type of systematic error, namely pulse-length error (PLE), 
occurs when the radiation (assumed to be a rectangular pulse), is resonant with the target 
spin transition but whose application-time or amplitude differs in an unknown (but fixed 
in time) quantity from the ideal value. Another type of systematic error is off-resonance 
error (ORE). Off-resonance errors arises when the frequency of the control radiation is 
unexpectedly not on resonance with the spin transition and spin dynamics proceeds with an 
unknown (but fixed), detuning parameter. 

In this paper, we describe schemes of generating an entangling gate between the elec- 
tronic and nuclear spins in the NV centre in diamond which are robust against PLE and 
ORE systematic errors. In this letter, we investigate effects of such errors in a sequential ap- 
plication of rectangular pulses of microwave and radio-frequency. We compare the fidelity of 
such pulses with their more robust counterparts: composite pulses, and numerically derived 
GRAPE pulsesjl2|. The sequential pulses have previously been used to generate two-qubit 
entanglement 7|, while composite pulses are known to be capable of correcting systematic 
errors in nuclear magnetic resonance (NMR) experiments 13|]. In the next section, we de- 
scribe the Nitrogen- Vacancy (NV) system and introduce the sequential pulses required to 
create a particular entangled state in this system. Section HI examines the decrease in gate 



|2) = |l)e |0)„f ♦ 



3A 



^02 



|0) = |0)e® |0), 



1 
1 

1 > 




\5r- 



























|3) = |l)e® |l)n 



|1) = |0)e® |l)n 



FIG. 1: Energy level diagram of the ground state of the NV centre in the absence of an external 
magnetic field. The |0) — |1), and |2) — |3), enery splitings are due to the hyperfine interaction 
between the NV electron and the nuclear spin of the 13 Carbon atom. The |0) - |2), |2) - |3) 
transitions are driven by MW, RF radiation observing selection rules. 

fidelity due to systematic errors when we adopt the sequential, robust composite and robust 
GRAPE pulse control. Results of the numerical simulations and performance comparisons 
between the composite and GRAPE pulses against PLE and ORE systematic errors are 
discussed in Section IV, together with conclusions. We find that GRAPE pulses are more 
robust than the sequential and composite pulses, and faster than the composite pulses. 



II. ENTANGLEMENT GENERATION 



To describe the NV, we follow a model described in detail in [lOj and concentrate our 
discussions on the fine and hyperfine structure of the ground state of the NV centre coupled 
to the ^^C atom. When there is no external magnetic field aligned with the quantisation 
axis of the NV centre, the system has four levels in the ^A manifold as depicted by Figure 
1. This is due to the degeneracy of the = ±1 electronic spins (5 = 1), of the centres 
and the interaction with the nuclear spin (/ = 1/2), of the Carbon atom. For clarity and 



consistency, we adopt the same notations for the four spin levels as those used in 
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The states {|0)e, |l)e} are eigenstates associated with the electronic spins {m^ = 0, = —1} 
while the states {| t)„, | |)„} correspond to the nuclear spins {m/ = — l/2,m/ = +1/2}. 
From [10|, the |0) — |2) and |2) — |3) transitions are found to be uJo2 ~ 2.88 GHz and 
tuos ~ 130 MHz respectively, while uqi ~ 2 MHz. In general, the two frequencies Um and Ur 
can be set slightly off resonance as parameterised by detunings 6m and 6r- These detunings 
have unknown but constant values in the case of ORE. 

The Hamiltonian of the system described in Figure 1 can be expressed as 

(fi™e^-™*(T2o + f^.e^-^*a23 + H.c.) , (5) 

where dpg = \p){q\, flmi^r) are the Rabi frequencies driving the MW (RF) transitions. The 
first line of (5) is the self-energy of the system relative to the ground state |0), while the 
reminder of the Hamiltonians describes the interaction of the external radiation with the 
system. We now move to the interaction picture (IP) defined by Uo{t) = exp{—iHot), with 
Hq = aa22 + ba^s + can + daoo,, to obtain Heff = UHU"^ — iU'^U ^ 

Heff = {uJo2 - a) 0-22 + {(^03 - b) 0-33 
+ (t^oi - c) 0-11 - d&oo 

{QmCr20 + ^r^23 + H.C.) , (6) 

where we have set a — 6 = Ur, a — d = Um- Defining / = aoo + <5"22 + o"33 + (5"ii, (T^° = 
—(Too + o"22, (^T = o"22 — 0-33, = o"33 + <5"ii,, givcs the effective Hamiltonian a new form 



expressed as 



1 
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Sr — COoi — a + C) (Tl 
^m- + + UJqi + a - c) Crf 

—SuJm — 26m + + SuJoi + 3a — 3c) 
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-- + H.C.) - - {QrCr23 + H.C.) . 

The terms having / and in the above expression vanish by taking a 
cjoi,, so that, 
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= 5. 



where express the control pulse Vt 

m,r — ^m,r Gxp (z0m,T-) and we have chosen 6^ 
We note that Um,r and 6m,r are a real, and describe the control amplitudes and control 
phases, and use a^'' = cXpg + a^p, a^^ = i {cTpg — CTqp)- We can clearly see from ([H]) that the 
Hamiltonian of the system is reduced effectively from the four-state system to a three-state 
system involving only the states |0), |2) and |3), which can be expressed as 
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We now consider the sequential application of MW and RF unitaries via rectangular control 
pulses. We wish to generate the entangling gate |\E'/) = -R23(7r)-Ro2(7r/2)|0) = JJrUm\^)- We 
consider the resonant case (5 = 0), and we first apply a MW pulse on the |0) ^ |2) transition 



for time tm = ^u'} to get f/„ = e^'^y „ and = f/^|0) = ^(|0) - |2)). Switching off 
the microwave radiation and subsequently applying a radio-frequency pulse with a phase of 
Or = |, for a duration of tr = ttu^^ on resonance with the |2) ^ |3) transition produces the 

^ TT ■ - 23 

gate Ur = e2*°"a , which transforms the state \iprn) into another superposition state \iPb), 



^(|0)-|3)) = -^(|0,0)-|1,1)). 



(10) 



The state is an entangled Bell state in the coupled system of electronic and nuclear 
spins. The sequential application of the on-resonance microwave and radio-frequency pulses 
generates a sequential unitary gate Usq, 
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which takes a total time t 
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III. SYSTEMATIC ERRORS 



To study the effects of systematic errors on the gate fidelity, we first quantify the PLE 
and ORE by two error fractions ef = (T' — T)/T, and eg = 6/ A, respectively, where —1 < 
{^f, Cg} < 1- We use T and T' to denote the ideal and non-ideal pulse application times and 
A is taken to be the fixed maximum amplitude of the microwave or radio-frequency pulses. 
In the presence of off-resonance error we consider the case when in the sequential pulse, 
maxj-Um} = max{Uf.} = A. 

When the sequential pulse suffering from pulse-length error is applied to the system the 
actual gate executed is highly dependent on the error fraction e/, and we can write. 



frf _ 77/77/ _ pn/2{l+ef)iaf n/i{l+ef)iaf 
'-^sq '-^r '-^m ^ ^ 



[12) 



We use the standard gate overlap fid elity F to measure the closeness between the generated 
gate, say f/a, and the target gate JJi [14 1, 
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One can show that the fidelity of the actual gate f// with respect to the ideal one tfsq is 
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and is numerically plotted in Figure 3. The quadratic term of e/ in ([T3j) significantly reduces 
the fidelity. Hence, It is very desirable to find other types of control pulse that can suppress 
and possibly eliminate the quadratic term of e/ from the gate fidelity. 



In the presence of the off-resonance error, the actual gate implemented by the sequential 
pulse is 

= U^e.me,), (15) 

where 

1/9 = (16) 

[/9 = e-t*^«(^'°+^'')+f'^'". (17) 

Due to the complexity of the analytical expression, the fidelity of the gate tJ^^ is calculated 
numerically and plotted in Figure 5. 

A. Composite pulses 

Originally invented in NMR by Levitt and Freeman [15], composite pulses offset the effect 
of systematic errors by replacing single quantum operations with several quantum operations 
which are designed to cancel out systematic errors. Composite puls es have been developed 
to implement robust arbitrary single qubit gate operations 13|, [l6| . In our case, we apply 
two particular types of composite pulses, (A) Broad Band Number 1 (BBl) composite pulses 
to correct for PLEs and (B) compensation for off-resonant errors with two CORPSE pulses. 
To apply these composite pulses, we first define generalised forms of the pulses subject to 
imperfect microwave and radio-frequency controls. In the presence of pulse-length error, 
these generalised forms are 

ULijm, Om) = e^C^o^e^^'^+^i'i^-^f (18) 
W/ {Tr,9r) = e*i('=°'^^'-°"^^+^™^'-'^a^)"'-(^~'=/)'^'-. (19) 

On the other hand, the generalised quantum operations in the case of off-resonance errors 
are 

W{Tr,9r) = e~*(^'^~5(™^^™°"^^+'^™^'"°'a^))"'-'^'-^ (21) 

where Z = af + of. 



We formulate the composite pulse analogue of the sequential gate ( fT2i) by replacing [/^ 
and jJl with their composite counterparts. For pulse-length errors we use BBl composite 
pulses to replace Uj^ and f// by C^^ and C/ 



Ci^ = Ul{7:|^,7^|2)Ul{'K^Mlx)Ui{2'KX^2^x^)^^^ (22) 
Cl = W/(7r/2,7r/2)W/(7r, 1.087r)W/(27r,2.247r)W/(7r, 1.087r)W/(7r/2,7r/2), (23) 

where and are expressed in the units of and u~^. One can verify that in the 
absence of PLE, we have C/(e/ = 0)C^(e/ = 0) = Usq. The composite pulses C-^ = C/C^^j, 
takes a total time r-^ = Qttm"^ when both real amplitudes Um = Ur are the same. The total 
time is exactly 6 times longer than that of the sequential pulse. The fidelity of the BBl 
composite gate, C-^, plotted in the Figure 3, is numerically calculated through the following 
equation, 



Tr ( UlqUsg 



(24) 



The range of good fidelity (defined for F > 0.9) has expanded. In contrast with the sequential 
gate, where the good fidelity holds approximately only for |ef| < 0.4, the BBl composite 
gate maintains the good fidelity for |e/| < 0.7. 

In the case of ORE, we replace the gates f/^ and IJ^, with these CORPSE counterparts 

H, 

= W^(2.14vr, |)W^(1.777r, -|)W^(0.14vr, |), (25) 
C9 = U^{7n/3, |)W,n57r/3, —p^T^/S, |), (26) 

to produce the complete gate = CfC^, which takes a total time of ~ 5.59 x tgg, nearly 
six times longer than the sequential pulses assuming Um = Ur- We plot the fidelity of the 
gate numerically calculated based on Equation (fT3ll in Figure 5, and find that the CORPSE 
pulses do not correct the off-resonance error. This is due to the Autler-Townes splitting of 

n 

the NMR transition due to the ESR excitation of the MW transition p/7| . Thus one cannot 
use the naive CORPSE pulse in each of the ESR and NMR sequential pulses to correct for 
ORE. 



B. Robust GRAPE Pulses 



We now explore whether one can obtain robust operations in a shorter time than the 
composite pulses by simultaneous irradiations both the MW and RF transitions. Rapid and 
robust control quantum control is important as the total gate duration significantly deter- 
mines the number of quantum operations that can be performed before decoherence degrades 
the quantum coherences. Initially proposed for NMR experiments, the GRAPE algorithm 
produce pulses that minimise the time required to implement a target unitary operator 
12 1 and some quantum algorithmic elements [l^. GRAPE pulses have been experimentally 



n 

demonstrated in a single qubit trapped ion system . A GRAPE based scheme is also pro- 
posed to control a coupled Josephson qubit system [l9|], and later extended to control open 



quantum systems in the Markovian domain 



20|. It has also been applied in implementing 



high-fidelity single qubit operations in a noisy environment due to random telegraph noise 



in superconducting solid-state qubits 



2l| 



We briefiy summarise how the GRAPE algorithm works 121] and how it can achieve time- 
optimal control and robustness against systematic errors. We start by writing the unitary 
evolution under the Hamiltonian (IHl) in the form of 



ij = T /e-<^=+^*-"^«^0*(it, 



(27) 



where Uk{t) and Hk are the control pulses and control Hamiltonians, expressed as follows 

■^cose^(t), "'^^^^ 



= cos e^it), u^it) 

Usit) = — ^^^^ COS 6'r(t), U4{t) 



2 sin6'^(t), 
-^^sin Orit), 



(28) 



and {Hi,H2,Hs,H4} = {a^", a^^, cr^ }, and Hg is a drift Hamiltonian, Hs = 
|5 (o"^° + (T^^) = ^6Z, and with the initial condition, U(t = 0) = /. We wish to numer- 
ically optimise the controls Uk{t) in a particular application time t = T, such that U{T) 
approaches a target gate Ut- This is equivalent to maximising a performance function P, 

2 



Tr[/|f/(T) 



(29) 



One considers discretising time At = T/N, where is a number of time steps/bins, and 
during each time bin the control u^s are constant. We can approximate fl27p . 



f/ = f/^f/^_i...f/i. 



(30) 



which gives an approximation of the performance function to be 

P = Tr {UlU nU N -I ■■ A) ^ 

Tr((f/i...t/;v-if/7v)^f/r). (31) 

From the gradient g{j) = 5P/Suk{j), to 0{At), is written as 

g{j) = ~2Ke{Tr{tAtA]HkB,)Tr{B]A,)), (32) 

where Aj = Ujj^_-^^...UlfUT and Bj = Uj...Ui, and the performance function P always increases 
if we update 

Mfc(j) -^Mfe(j) + e^(j), (33) 

where a small step size e is used. It is also necessary to add an additional term gmax{j), 

gmaxU) = -2apUk{j)At, (34) 

to the gradient g{j), to penalise excessive microwave or radio- frequency power. 

The pulse-length and off-resonance errors are incorporated into the algorithm by replacing 
the ideal Uj in Equation ( l30l) with Uj and f/J written as 

The performance function to be optimised is the average performance function over set of 
error fractions, and defined as 



Tr((f/;^-...?7;/l^f/;/-)tf/^), (37) 



for the case of PLE or ORE (denoted by superscripts f or g), where N{tf^g) is a number of 
elements in the set of error fractions ef^g we wish to optimize over. This modification allows 
us to optimise the control pulses Uk only for a certain range of error fractions, for example 
—0.2 < ef^g < 0.2. However, the real performance of the optimised pulse is checked through 
the gate fidelity as defined by Equation (fT3l) . 




Time (s) 

FIG. 2: The real amplitudes and phases of the microwave and radio-frequency pulses engineered 
by the GRAPE algorithm to create an entangling gate which is the robust against the pulse-length 
errors (PLE). The application time is 67r/xs and the maximum real amplitude is 1 MHz. 

IV. DISCUSSIONS AND CONCLUSIONS 

In our numerical simulations, the modified GRAPE algorithms explained in the previous 
section are executed using Matlab to produce a set of optimised microwave and radio- 
frequency pulses presented in Figure 2 and 4. Both sets of pulses have N — 400 and 
max{um} = max{ur} = 1 MHz and this limit is realistic in current NV ESR/NMR ODMR 
experiments. We found that by simultaneous irradiations we are able to reduce the time 
taken to achieve robust operation with respect to both the pulse-length and off-resonance 
errors. We manually re-optimise the pulses by adjusting the range of the error fractions and 
the duration several times until the gate reaches a considerably high fidelity. 

In both cases of pulse-length and off-resonance errors, the gate fidelity of the optimised 
pulses via GRAPE outperforms those of the sequential pulses and the composite pulses for 
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FIG. 3: Fidelity plots for the sequential pulse (Line 1), the BBl composite pulse (Line 2) and the 
GRAPE pulse (Line 3) against PLE. 

every value of the PLE and ORE fractions, as shown in Figure 3 and 5. However, the BBl 
composite pulse does better than the sequential pulse in the case of pulse-length error. It 
is interesting to note that the CORPSE composite pulse, which is robust against ORE in 
a single two-level system, is no longer robust against the ORE in our essentially three-level 
system. Hence, a further investigations would be needed in order to develop composite 
pulses that can tackle the off-resonance error in a three-level system. This is beyond the 
scope of this paper. 

While the maximum real amplitudes of the optimised pulses used in our analyses are not 



technologically difficult to implement |8|, l9|, lul] , the fast controlling apparatuses required to 



rapidly change the pulse amplitudes and phases is quite demanding since they should be 
able to operate on nanosecond timescales. 

In conclusions, we have numerically optimised the microwave and radio-frequency pulses 
required to create the entanglement in the system of single NV centre coupled to the nearest 
Carbon atom through the modified GRAPE algorithm. We have found that the entangling 
gate created by such optimised pulses is more robust against systematic errors and has faster 
implementation time than that required by the corresponding composite pulses. Reasonable 
extensions of our work would be considering the simultaneous presence of PLE and ORE in 
the system and taking into account decoherence processes. The latter one certainly needs 
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FIG. 4: The real amplitude and phases of microwave and radio-frequency pulses engineered by 
the GRAPE algorithm to create an entangling gate which is the robust against off-resonance errors 
(ORE). The application time and the maximum real amplitude are the same as those in the case 
of PLE. 



the system to be modelled by an open quantum system analysis through master equation. 

This work was supported under the European Commission FP6 1ST FET QIPC project 
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FIG. 5: Fidelity plots for the sequential pulse (Line 1), the CORPSE composite pulse (Line 
and the GRAPE pulse (Line 3) against ORE. 
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